Experimental evidence for a two-band superconducting state of NbSe2 single crystals 
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We report on measurements and a detailed analysis of the reversible magnetization of supercon- 
ducting NbSe2 single crystals. By comparing the experimental data with Ginzburg Landau theory 
we show that superconductivity in NbSe2 cannot be explained by an anisotropic single-band, but 
by a multi-band scenario. Applying a simple two-band model reveals the basic mixed-state pa- 
rameters, which are quite different in the two bands. We identify a strongly anisotropic band that 
determines the properties at high magnetic fields, and a second almost isotropic band that dominates 
at low fields. Our method is well suited for distinguishing anisotropic single-band from multi-band 
superconductivity in various materials. 

PACS numbers: 74.25.Ha,74.25.Op,74.70.Ad 



I. INTRODUCTION 

NbSe2 is certainly one of the most intensively stud- 
ied superconducting materials for several reasons. First, 
its transition temperature (Tc) is about 7K and its up- 
per critical field not much larger than 4 T perpendicular 
(B^2i c-direction) and 12 T parallel {B'^^i ab-direction) 
to the Nb planes."^ Accordingly, most part of the super- 
conducting phase diagram is accessible to experiment, in 
contrast to many other materials. Furthermore, large 
high quality single crystals with almost negligible vor- 
tex pinning effects can be grown. Introducing a small 
amount of disorder (e.g. by particle irradiation) may lead 
to the emergence of the well known fishtail effect,^ which 
is still discussed a lot by the superconductivity commu- 
nity. Moreover, NbSe2 was the first material, in which 
scanning tunneling microscopy was successfully employed 
for observing vortex cores or distributiona^i^ and it is still 
widely used for such investigations. Finally, the charge 
density wave state, formed below about 33 K,^ allows 
studying the effect of competing order parameters, which 
is an important issue for high temperature superconduc- 
tors. 

Recently, NbSe2 was suggested to be a two- or multi- 
band superconductor, which was confirmed by several ex- 
periments (e.g. Refs. l6l-[lll). but usually an alternative 
interpretation of the results in terms of an anisotropic 
s-wave single-band scenario could not be excluded. In 
this paper we provide further evidence for the two-band 
scenario. The field dependence of the reversible mag- 
netization - M{B) - is compared with Ginzburg Lan- 
dau theory, showing that M{B) cannot be reliably de- 
scribed by the anisotropic single-band scenario, but by 
multi-band superconductivity. Evaluating the curves re- 
veals the anisotropy and other superconducting parame- 
ters which are not only temperature but also significantly 
field dependent. 

More generally, following the discovery of superconduc- 
tivity in MgB2, multi-band scenarios were announced for 
a lot of materials including Fe-pnictides, borocarbides. 



heavy fermions, and even cuprates 



12-16 



etc. However, 



the experimental evidence was often based on poor argu- 



ments (like, e.g., the observation of a temperature depen- 
dent anisotropy) and the anisotropic single-band scenario 
could rarely be eliminated. Therefore, a method for dis- 
tinguishing between the two scenarios is clearly desirable, 
which is provided by our approach. At the same time, 
all basic mixed state parameters can be assessed in this 
way. 



II. EXPERIMENT AND EVALUATION 

Two NbSe2 single crystals grown by a standard chem- 
ical vapor transport technique with sizes a x b x c c^i 
2.25 X 2 X 0.1 mm'^ and 1 x 1 x 0.1 mm^ were investi- 
gated. Both have a Tc of 7.15 K and a small transition 
width of about O.IK. The smaller crystal was mainly 
used for confirming the results of the larger one. 

Measurements of the magnetic moment parallel to 
the applied field (m) as a function of the applied field 
{m^a. < 7T, ^0 = 47r X 10"'^ T m A"^) were carried out 
at temperatures from 2 K to about 10 K in our SQUIDs. 
The curve above T^ (at 10 K) was subtracted from the 
superconducting signal to get rid of the temperature in- 
dependent background signal (which was non-significant 
in most cases). In all measurements, most parts of the 
magnetization curves were reversible. Therefore, the re- 
versible magnetization could be obtained directly from 
M = m/V {V is the sample volume). In the case of irre- 
versibility, the reversible signal had to be calculated from 
the irreversible parts at increasing (m+) and decreasing 
{rri-) field via M — {m++m_)/2V . Minor corrections to 
the field coming from the critical current were considered 
numerically as described in Ref.llTl Finally, the magnetic 
induction B was evaluated via B = ij,o{Hs^ — DM + M), 
where D is the demagnetization factor of the sample, 
and M{B) obtained. Fitting M{B) from single-band 
Ginzburg Landau (GL) theory to the experiment results 
in two parameters, the upper critical field (Sc2), which is 
the field where M{B > 0) vanishes for the first time, and 
K, the GL parameter. For simplicity, the GL results were 
approximated by simpler (but quite accurate) equations 
given in Ref. |l^ (we used Eq. 23 of that reference leading 




FIG. 1: (Color online) The reversible magnetization from ex- 
periment and GL theory. The open circles illustrate the ex- 
perimental data of NbSe2 for i/a || c (panel a) and Ha. \\ ah 
(panel b) at 4.2 K. The dashed lines present the best single- 
band GL fits to the whole field interval, the solid lines to the 
low- or high-field region only. The inset of panel (a) shows a 
measured hysteresis loop from which the reversible data has 
been evaluated. The inset of panel (b) presents the reversible 
magnetization of a VaSi single crystal at 13.5 K (open sym- 
bols) and the best GL fit (solid line). 



to deviations from the exact behavior by less than 3.5%). 
We point out that other approximations or even numer- 
ical solutions of the reversible magnetization M{B) may 
be applied, which are obtained from GL theory or other 
models. 



III. RESULTS AND DISCUSSION 



Fitting the reversible magnetization by single- 
and two-band models 



Figure[T]presents results of the experiments and the fit- 
ting procedure. The full circles in the inset of panel (a) 
show the low-field part of the measured magnetization 
loop m/V for Tfa || c at 4.2 K. A significant hysteresis 
occurs only at very low fields (below about 20 mT, which 
is ~ 0.01i3c2) and this region was excluded from the re- 
versible curves presented by the open circles in the fig- 
ure. The dashed lines show the best single-band GL fit to 
the whole superconducting part (from about 0.01-Bc2 to 
Bc2) and therefore represent the results we expect from 
a single-band superconductor. Here, k was the only fit 
parameter while Bc2 was taken from the point, where the 
experimental curves reach zero. Enormous differences be- 
tween experiment and the best GL fits are evident (e.g. 
Fig. [1]) at both field directions and all temperatures (2 
- 6.6 K). These deviations become even more apparent 



when fitting the GL theory only to the high- or low-field 
region (in the latter case, also i?c2 is fitted) as illustrated 
by the solid lines in Fig. [TJ We conclude that NbSe2 ob- 
viously does not behave like an anisotropic single-band 
superconductor . 

The GL theory holds strictly only sufficiently close to 
the phase transition and T^. At low temperatures and 
far from Bc2 some deviations are expected, but usually 
GL theory is quite successful over the whole field and 
temperature range. Nevertheless, to clarify this point, 
we fitted the single-band GL theory to magnetization 
curves of Nb and VaSi, which are single-band supercon- 
ductors. In all cases the GL fit showed reliable agreement 
with experiment over the whole field range, as illustrated, 
e.g., in the inset of Fig. [Tfb) for VsSi. Fair agreement 
between GL theory and M{B) data was also found in 
the more anisotropic superconductorsi^r^ YBa2Cu408, 
Ndi.85Ceo.i5Cu04-5, and YBa2Cu307_5. 

In the next step we analyzed possible two-band ef- 
fects. In the anisotropic single-band scenario, proper- 
ties like the energy gap or the coupling strength vary 
within one band, while in the multi-band modeP^ super- 
conducting charge carriers exist in two (or more) electron 
bands at the Fermi surface characterized by two (or more) 
sets of parameters. Two-band behavior can be observed 
if two of the parameter sets are significantly different, 
e.g. for the upper critical field or the anisotropy, and if 
both bands contribute substantially to the overall prop- 
erties. Additionally, interband coupling (pairing of elec- 
trons from different bands) is expected, which leads, e.g., 
to one single transition temperature of both bands. Even 
quite small interband coupling strengths would strongly 
mask the superconducting transition of the band with 
the lower values of Tc and Bc2, which, however, does not 
mean that the properties follow the behavior of a conven- 
tional single-band superconductor.— While the temper- 
ature dependence of most quantities is often quite sim- 
ilar for anisotropic single-band and two-band materials, 
the field dependence is rather different. To analyze two- 
band effects in the experiment, we separated M(B) into 
a high- and a low-field region and fitted both regions in- 
dependently by single-band GL theory (see solid lines of 
Fig. p. 

In the high-field region, good agreement was achieved 
in the range from about 0.6 B^2 to B^2 for -^a || c and 
from 0.5 5^2 to B^2 &!■ Ha \\ ab at aU temperatures. This 
indicates that superconductivity of the second band is 
suppressed at those high-fields and only the single-band 
properties of the band with the larger upper critical field 
(indicated by the superscript 'a') survive. Some effects of 
interband coupling cannot be excluded, but are expected 
to be small. Accordingly, the high-field GL fit is assumed 
to roughly provide the single-band parameters of the a- 
band (see Fig. [2 and Tab.I|. 

Turning to lower fields in Fig. [Jl the discrepancy be- 
tween the experimental data and the high-field single- 
band GL fit increases, which we ascribe to the influence 
of a second band with a smaller Bc2 (superscript 13). 



With decreasing field the /?-band becomes more domi- 
nant (particularly for if a II ab), but the contribution of 
M°'{B) does not become negligible. Therefore, single- 
band behavior is not expected and the low-field GL fit 
only provides effective parameters (index 'If') of the over- 
all superconducting behavior, although more indicative 
of the /3-band. The fit interval ranged from 0.05i?^2 to 
0.3-6^2 fo'" -^a II c, which was the largest interval that 
enabled good matching at all temperatures. The fit pa- 
rameters are more affected by reducing the interval than 
in the high-field regime, but the major findings are not 
concerned. To be consistent, we applied the same field 
interval in absolute values for H^, \\ ab, which roughly 
corresponds to Q.QIQB^^ to O.IB^^^ 

To come closer to the single-band properties of the 
/3-band, we subtracted the a-band GL fit from the ex- 
perimental data. However, pure single-band behavior of 
the remaining M{B) {M^{B)) is still not expected and 
indeed not observed due to interband coupling (which 
might dominate, where M^{B) approaches zero, i.e. near 
Bc2 of the /3-band) and possible additional superconduct- 
ing bands. Accordingly, the field interval, where single- 
band GL theory matched well was similar as for the effec- 
tive 'If regime before subtracting the a-band. In contrast 
to the 'If regime, however, the thermodynamic critical 
field {Be, see Tab. |l| of the /3-band is almost indepen- 
dent of the field orientation as expected for a single band 
(small deviations are ascribed to the interband coupling 
and experimental uncertainties). 



B. Mixed-state properties and anisotropy 

The results of the fitting procedure are illustrated 
in Fig. [2] and Tab. H B^f (a-band) follows well the 
behavior of simple conventional superconductors (full 
circles in Fig. [2la)). For extrapolating the data to 
OK the model Bc2{T) = Bc2(0)[l - (T/Tc)"]'' (lines in 
Fig. IHa)) was applied which matches best for a ~ 1.3, 
6-1, and S;;f (OK) ~ 4.1 T. The same model leads to 
B^^'^iOK) ~ 12.3T (i/a II ab). In the latter case, we 
observe a positive curvature of i?c2 near Tc, while -8^2" 
is linear in this regime. The corresponding GL parame- 
ters (k. Fig. [21b)) decrease strongly with increasing tem- 
perature. Extrapolation gives roughly k'^'°'{OK) ~ 30 
and K^-^iGK) ~ 17 as well as /{^"^^"(OK) ~ 86 and 

^ab,a(-ggj^) -30. 

The /3-band upper critical fields are also well fitted 
by the above model leading to much smaller results at 



1.6 T and B^^'^^iOK) 



1.9T. 



OK, namely B^^:l'{OK) 
Also the GL parameters decrease and show a less pro- 
nounced temperature dependence, i.e., k'^'^(OK) ~ 15 
and k'^'^(6K) ~ 15, ^^'^''^(OK) ~ 21, k^'^''^(6K) ~ 19. 
The effective parameters (which address the sum of the a- 
and /3-band) of the low-field regime deviate only slightly 
from the /3-band data (see Fig. [2] and Tab.[T|. 

The above results lead to the anisotropies 7bc2 = 




B,t/i?c^2 and 7k 



^ab 



/k" illustrated in panel (c) of 



FIG. 2: (Color online) The superconducting parameters of 
NbSe2 obtained from evaluating the reversible magnetization 
curves. Panel (a) displays the upper critical fields (-Bc2), panel 
(b) the Ginzburg Landau parameters (k), and panel (c) the 
anisotropies. The circles (with solid hues) present a-band, 
the squares (with sohd lines) /?-band, and the triangles (with 
dashed lines) effective low-field regime results. In panel (a) 
and (b) the open symbols refer to H^ \\ ab and the full symbols 
to _ffa II c, in panel (c) to the Bc2 (open) and the k (full) 
anisotropies. The lines are guides for the eyes. 



Fig. [21 At high fields (a-band) and temperatures be- 
low ~ 5K, 7bc2 — 3 and almost constant (open circles), 
i.e. close to the value usually reported for NbSe2,^ but 
drops to 2.3 at 7K. Such a temperature dependence cor- 
responds to the positive curvature of B^2iT) near Tc 
and is, in general, a characteristic of both, two-band— 
and anisotropic single-band^^ superconductivity. The 
anisotropy of k (full circles) is almost identical to 7bc2 
(open circles) in the high-field regime. At lower fields the 
anisotropy of the GL properties is considerably reduced. 
We obtain 7*^ ~ 1.3 (7302 ^ 1-2, 7« ~ 1.4) and 7" ~ 1.35 
(7Bc2 — 1-2, 7k — 1.5), thus almost isotropic behavior. 
The small differences in 7bc2 and 7„ could result from 
the larger experimental uncertainties in this regime. 

In summary, the anisotropy of NbSe2 changes strongly 
with the applied field from a large value (~ 3) at high 
fields (> 0.5i?c2) to almost isotropic behavior at low fields 
(< 0.3i?^2; ^ 0-l^c2 )■ It seems impossible to explain this 
behavior by a single-band model, i.e. a two- or multi- 
band model is needed. Our results suggest the existence 
of a strongly anisotropic (a) band with a large i?c2, and 
a second almost isotropic {/3) band with a lower Bc2 (or, 
to be consistent with two-band theory: both bands may 
have the same Bc2 due to interband coupling, but the 
properties of the /3-band are strongly suppressed at high 
fields). We further point out, that the anisotropies of 
Bc2 and k are equal or quite similar, when acquired at 
the same field and temperature. Accordingly, this holds 
also for the anisotropy of all further parameters, such 



TABLE I: Summary of mixed-state parameters of single crys- 
talline NbSe2 at OK in the a- and /3-band and the effective 
low-field ('If') regime. Minor inconsistencies with the GL rela- 
tions occur from extrapolating the experimental data to OK. 
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Bc2 (T): 


4.1 


1.6 


2.0 


BS (T): 


12.3 


1.9 


2.6 


Bk (mT): 


9.0 


11.4 


19.1 


Bet (mT): 


3.9 


7.7 


10.3 


B^, (mT): 


98 


67 


112 


Bt (mT): 


98 


64 


86 


Xab (nm): 


265 


215 


165 


Ac (nm): 


795 


355 


350 


S,ab (nm): 


9 


14 


13 


Cc (nm): 


3 


12 


10 


k": 


30 


15 


13 


^"^ 


86 


21 


21 


T- 


3.0 


1.3 


1.4 











as the magnetic penetration depth (A) or the coherence 
lengths (^). These anisotropies were often claimed to be 
quite different in the two-band superconductor MgB2. 
But detailed evaluations showed that the situation is ac- 
tually the same as in NbSe2, i.e. all parameters have the 
same anisotropy*2^ The misinterpretation in MgB2 usu- 
ally originated from measuring the anisotropies of differ- 
ent properties at different applied fields, e.g. that of £, 
at high fields (from the Bc2 anisotropy) and that of A at 
low fields. 

Further parameters were calculated by applying the 
well known GL relations (e.g. ReL \2dt ). namely ^, A, 
Be, and Bd (the lower critical fields). The temperature 
dependence of these quantities was fitted by the same 
(general) model as Bc2 and the results at OK are listed 
in Tab. ID Additionally, B^ of the whole sample was 
obtained from integrating the reversible magnetization 
curve^^ i.e. —fJ-ojQ ^ MdH = i?^/(2/xo), resulting in 
about 120 mT at K. 



C. Discussion 

Our results are consistent with most findings from lit- 
erature. For instance, the Fermi level of NbSe2 is known 
to be crossed by three electron bands, one is related 
to Se-4p orbitals and shows a rather three dimensional 
pancake like sheet, and the other two to Nb-4d orbitals 
each leading to two rather two-dimensional (cylindrical) 
Fermi sheets;^ There is general consensus that a su- 
perconducting gap opens in both Nb bands. Investiga- 
tions by ARPES'' at 5.3 K resulted in gap values of 0.9 - 
1.0 meV, but no gap was detected on the Se bands. On 
the other hand, scanning tunneling spectroscopy^^ re- 
vealed a broader range of gap values from about 1.4 to 0.7 
(or even 0.4) nieV close to OK. The band, on which the 
smaller gaps exist, was not clearly identified, but it was 
argued, that the smaller gaps may open on the Se bands, 
which would not necessarily contradict the ARPES mea- 
surements, since the smaller gap could be suppressed near 
Tc in a two-band model. ^ It is plausible to associate the 
a-band, which has the higher i?c2, with the larger gap 
(Aq) and thus the /3-band with the smaller gap (A/j). 



The fact that the superconducting properties of the p- 
band are nearly isotropic may favor assigning A^ to the 
more three dimensional Se-band. On the other hand, 
ab-initio calculations revealed that the electron density 
of states of the Se band at the Fermi level - z(0) - is 
only about 5 % of the overall valuej^i which makes a 
substantial contribution of the Se bands to the super- 
conducting behavior questionable. Moreover, applying 
the BCS relations for the condensation energy density 
E^ = -z(0)AV2 = -B2/(2/xo) with B^ from Tab. Ill 
results in similar z(0) values for both bands, as was ac- 
tually found for the two Nb bands in Ref. [27| and thus 
favors a scenario, in which the two Nb bands are respon- 
sible for the two-band effects. It may be concluded that 
some microscopic aspects remain unclear and need fur- 
ther investigations. We note further that our results are 
in fair agreement with those from a specific-heat study.— 
Fitting the field dependence of the specific heat data by a 
simple two-band model led to similar upper critical fields 
and anisotropies of the two bands as in our study. 

Finally, we wish to mention the remarkable qualita- 
tive resemblance of the two-band properties in NbSe2 and 
MgB2, for which a similar analysis^^ also demonstrated 
a highly anisotropic band (cr-band) being dominant at 
high magnetic fields and a quite isotropic low-field band 
(TT-band), with similar values for the upper critical fields 
and the anisotropies, but a larger transition temperature 
(of almost 40 K). We point out that the mixed-state pa- 
rameters obtained by our approach in MgB2 are in good 
agreement with results derived from other methods, even 
for the band with the lower Bc2- 

There is also cumulative evidence for multi-band su- 
perconductivity in the faraily of the Fe-pnictide super- 
conductors (e.g. Refs. [l3ll3l ). Contrary to NbSe2 and 
MgB2, however, the anisotropy of the upper critical field 
(i.e., the high-field value of 7) seems to be smaller than 
that at low field a^^'^° (i.e. the anisotropy of the pene- 
tration depth or of the lower critical field). This would 
indicate that the band with the larger upper critical field 
is less anisotropic. More recent results suggest, however, 
that a proper description of experimental data is only 
possible by taking three or more bands into account i^i^^ 
Further insight might be gained when reversible magneti- 
zation curves over the whole field range become available 
for these materials. 



IV. SUMMARY 

In summary, we have shown that the reversible mag- 
netization of NbSe2 cannot be reliably explained by a 
standard anisotropic single-band scenario. Applying a 
simple two-band model reveals the anisotropy and other 
mixed-state parameters as a function of magnetic field 
and temperature. We found a strongly anisotropic band 
(7 ~ 3) that dominates at high fields and a second almost 
isotropic band relevant only at low fields. The method 
can be applied to other materials to distinguish between 



single-band and multi-band superconductivity and to ac- 
quire at the same time all basic mixed-state parameters. 
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